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Abstract 

The existence of light hidden sectors is an exciting possibiHty that may be tested 
in the near future. If DM is allowed to decay into such a hidden sector through 
GUT suppressed operators, it can accommodate the recent cosmic ray observations 
without over-producing antiprotons or interfering with the attractive features of the 
thermal WIMP. Models of this kind are simple to construct, generic and evade all 
astrophysical bounds. We provide tools for constructing such models and present 
several distinct examples. The light hidden spectrum and DM couplings can be probed 
in the near future, by measuring astrophysical photon and neutrino fluxes. These 
indirect signatures are complimentary to the direct production signals, such as lepton 
jets, predicted by these models. 
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1. Introduction 



The existence of a low energy liidden sector, weakly coupled to the Standard Model (SM), 
is an exciting possibility that will be tested by upcoming experiments. Hidden sector parti- 
cles can be produced in high energy colliders, as stressed in the context of 'Hidden Valley' 
models [1] and models where gauge kinetic mixing results in 'lepton jets' |2ll3]. Such hidden 
sectors can also be probed with low energy e^e^ colliders and fixed target experiments [i4j. 
Here, we point out that the existence of a low energy hidden sector, together with weakly 
interacting DM (WIMP) and gauge coupling unification, implies the generic possibility that 
DM may decay directly into the hidden sector through operators suppressed by the GUT 
scale. These decays, followed by decays into SM particles through kinetic mixing, provide 
the intriguing possibility of using astrophysical observations to study the hidden sector spec- 
trum, complementing direct production experiments. This decaying DM framework provides 
a simple and natural explanation for the recent cosmic ray (CR) anomalies [3], while avoiding 
the tensions and pitfalls of many previously proposed models. 

A DM explanation of the electronic CR excess requires a DM mass greater than a 
TeV [6l [7], and predominantly leptonic production [8J. Consequently, the vanilla MSSM 
WIMP scenario is disfavored, and many new models have been proposed, bifurcating into 
annihilating models P and decaying models [H], [121 lEl CH [E] . Annihilating models are 
difficult to reconcile with the FERMI and HESS CR data, because the softening of the spec- 
trum above a few TeV requires an annihilation cross-section (9(1000) times larger than that 
of the standard thermal WIMP [SIE]- Such a large cross-section is in tension with constraints 
from photon and neutrino measurements from the Galactic Center (GC) [SI [3, HEl [17] , extra- 
galactic emissions [18] , and the CMB [19] . There is also model building tension for achieving 
such a large cross-section. Possible mechanisms include non-perturbative Sommerfeld en- 
hancements [201 [H], or a resonance [22] • In the latter case, a very narrow resonance 
and degenerate states are required, while in the former, either large (> 1) gauge or Yukawa 
couplings to the light mediator or tuned parameters are necessary [9J. As we discuss below, 
the required large couplings conflict with a need for Yukawa interactions that generate a DM 
splitting, necessary in many models to avoid constraints from direct detection [23]. Indeed, 
the mechanism that generates the splitting typically opens up new annihilation channels that 
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can parametrically dominate at freeze out. As a consequence, in order to achieve the correct 
rehc abundance, the couphngs responsible for the Sommerfeld enhancement are constrained 
and cannot produce a large enough enhancement. 

Decaying models replace the need for a large annihilation cross-section. Since the DM 
lifetime is much longer than the age of the Universe, its decays do not affect the attractive 
features of the thermal WIMP and leave no signature on the CMB radiation. Moreover, 
constraints from the GC or subhalos are easily evaded [12], since the emission rate depends 
on one power of the DM density, p, as opposed to the dependence in the annihilating 
case. Interestingly, the correct lifetime to explain the anomalies, (9(10^^ sec), is obtained if 
the decays are induced by dimension-6 operators suppressed by the GUT scale [TT]. Still, 
it is non-trivial to construct a decaying DM model that does not over-produce antiprotons, 
and many existing models are fine tuned or have small and ad hoc parameters. 

In this paper we study a new and natural class of models, where DM decays into a light 
hidden 'dark sector', with gauge group Gd- Working in the supersymmetric framework ap- 
propriate in the context of GUTs, the dark sector has a stable mass gap at the GeV scale, and 
communicates with the supersymmetric SM (SSM) through kinetic mixing [2l|. The GeV 
gauge bosons decay into light SM fermions, explaining the lack of antiproton production 
The dark sector is close in spirit to the models discussed in [HI [2] • Nonetheless, it is more 
general in the sense that the DM may or may not be charged under Gd and/or the SM. This 
opens the door for a wider range of models and is potentially simpler. Dimension-6 decay 
operators appear naturally, and are expected to be present at low energy unless forbidden 
by global symmetries. For related work where DM decays into light states, see [I^. 

Models of the type studied here involve several scales. Physics at the GUT scale, Mqut, 
is responsible for producing the decay operators. More formally, in the limit Mqut ~^ oo, 
the DM is completely stable due to a preserved global symmetry. Fields at the GUT scale 
then break that symmetry, inducing the required decays. It is important that dimension-5 
operators which would trigger a fast DM decay are not generated. Below we show sev- 
eral mechanisms that prevent such operators from showing up at low energy. The TeV 
scale generates the DM mass which can be naturally related to the supersymmetry (SUSY) 
breaking scale, thereby avoiding the usual /i-problem. The GeV scale which controls the 
branching fractions of the DM decays into SM fields, is generated either by communicating 
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supersymmetry breaking to the dark sector indirectly through the SM ^ or through D-term 
mixing [31 EZj- Finally, sphttings between DM states may be required to avoid direct 
detection. Such splittings are naturally of order an MeV, thereby accommodating the in- 
elastic DM (iDM) [28j and eXciting DM (XDM) scenarios j29j . Below we study mechanisms 
that can appear at each of these scales, stressing the modular nature of such models, which 
significantly simplifies the model building. 

The models studied here predict distinctive signatures in many upcoming experiments, 
and unique indirect signals which will complement the direct production experiments men- 
tioned above. For instance, if the dark sector is approximately supersymmetric, or if the 
dark gaugino is lighter than the dark gauge boson, m;^^ < m^^, it typically decays into a 
gravitino and a SM photon. Such primary photons will show up as sharp features in the 
measured flux. If the DM is also charged under the SM, its decays are accompanied by 
primary neutrinos, again admitting a sharp and hard spectral feature. In the corresponding 
annihilating models, these decay channels are excluded due to the excess of primary photons 
or neutrinos produced, for example, at the GC. Both possibilities are studied in [30], where 
it was shown that current and future experiments will have the ability to measure these sig- 
natures and thereby differentiate between the annihilating and decaying DM scenarios. In 
sections |2] and [3] we provide detailed examples that illustrate the presence of these signatures. 

The paper is organized as follows. In section |2] we discuss the tools for constructing 
decaying DM models. We first list the dangerous pitfalls of these models, and then discuss 
solutions, organized by energy scale. In section [3] we apply these tools to study four distinct 



example models. In 3.1 we show the simplest U{1) model, which is UV completed in 3.2 In 
3.3| we construct a model where the DM is charged under the SM and decays into primary 
neutrinos, and in section 3.4| we demonstrate how one can evade direct detection without 
splitting the DM multiplets. In section |4] we discuss the cosmology of these models. In 
particular, we show that a supersymmetric dark sector can have long lived gauginos which 
decay into photons, without violating constraints from big bang nucleosynthesis (BBN). We 
conclude in section [5j In appendix |X] we revisit the symmetries of the four models, showing 
that these forbid the presence of any dangerous operators. 



4 



2. Tools for Modeling Decaying Dark Matter 

In this section we describe our strategy for building models of hidden sector decaying DM. 
After briefly introducing our framework and notations, we list several potential dangers for 
models of this type, which arise from cosmological and experimental constraints. We then 
introduce a series of model building tools, organized by energy scale, that address these 
dangers and can be used to build viable models. We stress that these tools are modular, 
and can be used to construct a variety of models. We demonstrate the use of these tools to 
build some example models in section |3| 



2.1. Framework 

We consider models where weak-scale DM, x, decays into a hidden sector with a gauge group, 
Gd, through a dimension-6 operator suppressed by the GUT scale, Mqut- We take this 'dark 
sector' to be weakly coupled with a GeV mass gap, in resemblance to the annihilating models 
proposed in Ref . |9] . Throughout this paper we work in the supersymmetric framework which 
comes naturally with GUT models, and can stabilize the GeV scale. Furthermore, we assume 
the breaking of supersymmetry to be mediated through gauge interactions, allowing for a low 
scale of mediation. This assumption can be somewhat relaxed, if the breaking is sequestered 
from the dark sector [31]. The dark sector consists of massive gauge bosons, 7^, gauginos. 
Yd, Higgses, hi, and Higgsinos, h^. We couple it to the SM through gauge kinetic mixing, and 
consequently dark sector particles decay through the mixing to SM particles. Due to the low 
dark gauge boson mass, it decays predominantly into light leptons. Within this framework, 
DM decays can naturally explain the PAMELA and FERMI measurements. Our notations 
are summarized in Fig. [ij 



2.2. Model Building Dangers 

• Dimension-5 DM Decay 

As we discuss in section |2l3, dimension-6 decay operators suppressed by the GUT scale 



induce DM decays with a lifetime of tq ~ 10^® sec, the correct timescale to account for 
the PAMELA and FERMI signals. Alternatively, dimension-5 operators suppressed by 
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GUT 



X, Y, {H) Decay operator 



TeV 



Xi, N„ (S) DM mass 



GeV 



7^, n, (h,) Gd broken 



Figure 1: We summarize our notations, organized by energy scale. X and Y denote GUT 
scale fields that are integrated out to generate dimension-6 operators that induce DM decays. 
We use {H) to denote a GUT scale VEV, which can partially break the dark gauge group, 



G'^ Gd, as demonstrated in section 3.2 The DM may be composed of multiple species, Xi, 
with mass at the TeV scale. This scale is naturally generated through the VEV of a singlet, 
S, that communicates with the SUSY breaking sector. Here, Ni denote electroweak scale 
fields that participate in the mechanism that generates a DM mass splitting. Such splittings 



can help evade the bounds from direct detection, as we discuss in section 2.6 The dark 
gauge group, with gauge bosons 7^, is entirely broken at the GeV scale by the VEVs of light 
Higgses. DM decays by dimension-6 operators into these GeV scale states. We use hi to 
denote light fields charged under the dark gauge group, at least some of which will receive 
VEVs, and we use n to denote a light singlet. 
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the GUT scale correspond to a lifetime of rs ~ 1 sec, and must be avoided. 
Sommerfeld Enhancement 

If DM is directly charged under the light dark sector, the annihilation cross-section 
is Sommerfeld-enhanced [9]. It is important that this enhancement is not too large, 
since there are various strong constraints on the annihilation rate. These include 
constraints from gamma rays and neutrinos from the Galactic Center and Galactic 
Ridge (GR) [HI El [17], diffuse gammas from extragalactic DM annihilations [18], and 
modified CMB radiation from DM annihilation during recombination 



Direct Detection 

There are strong limits from direct detection on models in which a weak-scale DM 
couples elastically to a light gauge boson that kinetically mixes with the photon. One 
finds a DM-nucleon cross-section of the order [9]: 

a,^^10-'^c^^(^Y(^)(^)-\ (2.1) 



,10-3; VO.Ol/ VI GeV/ 

where e parametrizes the size of the kinetic mixing. Current measurements rule out 
a cross-section of this size by 6 orders of magnitude [32]. There are also strong limits 
from direct detection on DM that couples elastically to the Z. For example, models 
where DM is the neutral component of an SU{2)w doublet are excluded by 2-3 orders 
of magnitude 



Inelastic Capture in the Sun 



As we discuss in section 2.6, one way to avoid the above constraints from direct de- 
tection is to split the mass between the DM states, 5Mdm ^ 100 keV, and couple 
inelastically to 7^ or Z [23] . It was recently demonstrated that if SMdm — 100 — 500 
keV, there are strong constraints on the inelastic capture of DM in the sun which 
is followed by annihilations into W'^W~ , ZZ, t^t^, or tt [31]. This constraint is 
particularly important if DM is charged under SU{2)w- 

Long-Lived GeV Scale Fields 

The dark sector may contain light long-lived fields, and one must make sure that their 
cosmology is safe. On the one hand, stable particles must not overdose the universe. 
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Qrch'^ < 0.1. On the other hand, the dark sector may contain long-hved particles that 
decay electromagnetically through the kinetic mixing. For such decays, lifetimes of 
order r ~ 10'' — 10^^ sec are constrained by Big Bang Nucleosynthesis [3S] and decays 
after recombination, r > 10^^ sec, are constrained by diffuse gamma rays |36j . 

• Long-Lived Colored Particles 

If the DM is charged under the GUT gauge group, then there is a colored component 
Xs- There are strong constraints on colored particles with lifetimes r > 10^^ sec as 
they form exotic atoms [37]. Xs niust therefore have a much shorter lifetime than x- 



2.3. GUT Scale: Decay Operators 

We consider models where weak-scale DM decays through dimension-6 operators suppressed 
by the GUT scale, into the dark sector. The GeV-scale dark fields then decay through gauge 
kinetic mixing to leptons. We focus on two possible scenarios, both of which include multiple, 
non-degenerate DM states: (i) One of the TeV fields receives a VEV, breaking part of Gd at 
the weak scale, and (ii) none of the states obtain VEVs and Gd is fully broken at the GeV 
scale. For scenario (ii), transitions between the TeV fields can be induced by the three body 
decay operators, 

I d'e x\x2h% , I d'e , I d'e xB/Wj. (2.2) 

GUT GUT GUT 

For the first two operators, xi ^i-nd X2 are both weak-scale with m^^ > ""^xi- Consequently, 
the DM is dominantly composed of X2 which generically has a larger density than xi- For 
the third operator x is a 5 of SU{5)sm- We will consider examples that generate each of 
these operators in section [3j The decay rate of these operators is given parametrically by: 

M^u ^-^,lo-.ec(ii^V7^%?^V. (2.3) 



167r2M^uT/ Vl TeV J \5 x GeV, 

This is the correct timescale to account for the PAMELA and FERMI signals, as was first 
noticed by Ref. [TT]. For scenario (i), two body decays will typically dominate. An example 



operator that we will consider in section 3.4 follows from inserting a (xi) VEV into the 



second operator of Eq. (2.2) 



As mentioned in the introduction, in the Mqut — * oo hmit, the DM is completely stable. 
This is typically achieved by a discrete symmetry under which Xi and Xi are charged. The 
superpotential at the GUT scale breaks this symmetry, destabilizing the DM. We demon- 
strate the existence of these symmetries in the models of section [3j Still, such symmetries do 
not ensure that the DM is sufficiently long-lived. Indeed, when integrating out GUT fields 
to generate the above dimension-6 decays, it is important to make sure that no dimension-5 
decay operators are generated. This can follow from symmetries at the GUT scale. For each 
specific model of section [3] we identify these symmetries in appendix |A| To demonstrate that 
this is possible, we now discuss two general mechanisms for generating dimension-6 decays 
that do not generate dimension-5 decays. One simple possibility is that the hidden sector 
gauge group is broken at the GUT scale, G'^ —>■ Gd, without breaking supersymmetry. By 
going to Unitary gauge and integrating out the massive G'^/ Gd vector superfields, it is sim- 



ple to check that dimension-6 decay operators, of the form of the first operator in Eq. (2.2), 
are generated in the Kahler potential ^38j. Moreover, if the DM and light Higgses have a 
canonical Kahler potential at the GUT scale, no dimension-5 terms are generated. We will 



discuss this in more detail for a specific example in section 3^ 

A second way to generate dimension-6 operators without generating dimension-5 ones is 
by coupling canonical GUT-scale fields to the DM in a chiral manner, 

W D MgvtXX + Xxh . (2.5) 

Integrating out X and X, and allowing for a weak-scale VEV for the DM, results in the 



dimension-6 Kahler potential operator in Eq. (2.4). It is straightforward to see from the 
equations of motion that no dimension-5 operators are generated in the superpotential or 
Kahler potential. More generally, global symmetries prevent quantum corrections from gen- 
erating dimension-5 decays in the Kahler potential, as we discuss in appendix [A} 



When the DM is charged under 5'[/(5)sm, as for the third operator of Eq. (2.2), one must 
ensure that its colored partner decays on a timescale shorter than the current age of the 
universe. For example, suppose that the DM is the neutral component of the doublet of a 
5-1-5. The model is viable if the triplet can decay through a dimension-5 operator that does 
not induce DM decays. This is straightforward to achieve since the triplet is typically heavier 
than the doublet at the weak scale, due to the RG evolution of their masses. Example triplet 
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Figure 2: A sample DM 3-body decay induced by one of the two last operators of Eq. (2.2). 
The DM decays to a GeV-scale gauge boson, gaugino, and a neutrino or the lighter field, Xi- 
The gauge boson decays through the kinetic mixing to a lepton pair and the gaugino decays 
through the kinetic mixing to a photon and gravitino, assuming that the gaugino is lighter 
than, or degenerate with, the dark photon. The resulting leptons can explain the PAMELA 
and FERMI excesses while the gamma rays and neutrinos lead to hard and sharp spectral 
features that can be probed by upcoming experiments [30] . 

decay operators include: 

^ ^ d'e , -J— [ d'e xioh , / d'e xBU , (2.6) 



where in the first operator the triplet partner decays into the DM while in the other two the 
triplet decays into a singlet, s, with m^^ < nis < tti^^. 

2.4- Weak Scale: Dark Matter Mass and Communicating SUSY Breaking 



As we discuss in the sections 2.5 and [5} sharp spectral features in the photon flux may exist. 



depending on the light dark spectrum. As a consequence, the low lying excitations, and 
indirectly the SUSY-breaking effects in that sector, may be probed in the near future [30] . 
Below, we briefly discuss the possible effects which may influence the spectrum. 

In our framework, the DM has a weak-scale mass. A GUT-scale one can be avoided by 
imposing a PQ or R symmetry that is spontaneously broken at the weak scale by a neutral 
scalar, S. The DM mass term then takes the form, 

y{S)xx- (2.7) 
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This is similar to the well-known /i-problem, and we present no new solution. Instead, we 
simply assume the coupling above, with a VEV induced by the SUSY-breaking sector. In 
principle, S may have a soft mass which arises from coupling to the SUSY-breaking sector. 
We distinguish between two cases, 

• X is not charged under Gd and couples to the light sector only through GUT 



suppressed couplings. Examples of such a scenario are given in sections |3.3| and |3.4 
In this case, SUSY-breaking effects are primarily communicated to the light sector 
through the kinetic mixing, as is worked out in [39j. The leading contribution to the 
soft mass squared of the light Higgses is generated as a threshold effect at the gauge 
messenger scale and is proportional to e^. 



,.nl . |m| ^ (100 MeV)^ (^) (|) 

Here gy is the hypercharge gauge coupling and is the soft mass of the right-handed 
selectron. As we show in the next subsection, this is parametrically smaller by one 
power of e compared to the supersymmetric mass squared of the dark vector boson. 
The corresponding contribution to the gaugino soft masses is even smaller [39] and 
may be neglected. The GeV scale, which we discuss below, is therefore approximately 
supersymmetric . 

X is charged under and may couple directly to the light Higgses. Such examples 



are given in sections 3.1 and 3.2 Here the SUSY-breaking effects can be communicated 
either through S or through the kinetic mixing as discussed above. Below, for simplicity 
we assume the latter. We stress that S can be naturally supersymmetric and still 
solve the /i-problem. This can be achieved for example through retrofitting [10]. We 
postpone the details of such a scenario to future work. If, on the other hand, S is 
accompanied by a soft mass, SUSY breaking effects in the light sector are expected to 
be of order GeV, and therefore dominate over the kinetic mixing contributions. 

In the above discussion we assumed the absence of TeV-scale messengers that couple 
to both the dark sector and SM. If such states exist, SUSY breaking is mediated as in 
gauge mediation and the supermultiplets are split at the GeV scale [2] . Finally, we note that 
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when the DM is approximately supersymmetric, both the fermion and boson components are 
cosmologically long-hved and constitute order one fractions of the DM rehc density. On the 
other hand, when there is large splitting within the DM supermultiplet, either the fermion 
or scalar component dominates the relic density, in a model-dependent fashion. The analysis 
that follows does not depend on the spin of the DM. 



2.5. GeV Scale: Breaking the Dark Sector 

In correspondence to the discussion above, there are two ways to naturally generate the 
GeV scale in the dark sector. One is with the use of D-term mixing which results from 
supersymmetric kinetic mixing |?T]. In such a case, the light dark sector is approximately 
supersymmetric, at the GeV scale. We review this mechanism below. The second way to 
generate the GeV scale is by communicating weak scale SUSY breaking as mentioned above. 
For simplicity, below we only consider a U{l)d model with the D-term mixing mechanism. 
The approximate supersymmetry in the light dark sector simplifies the analysis, since we 
do not need to consider GeV-scale soft terms. Nevertheless, we stress that this is only 
a simplifying assumption, which can easily be relaxed. Indeed, introducing GeV SUSY- 
breaking may change the low energy spectrum and consequently the astrophysical signatures, 
but does not affect the discussions below in a significant way. 

The GeV scale of our theory resembles that of [9l |2]. We assume that the SM and dark 
sector interact with each other through gauge kinetic mixing. The kinetic mixing between 
U{l)d and hypercharge is given by: 

''(feWdWy. (2.9) 

e is naturally of order 10^'^ — 10~^ and arises from integrating out heavy fields charged under 
both sectors. Supersymmetric kinetic mixing of this size automatically generates the GeV 



e 

~ 2 

\-3 in-4 



scale in the dark sector pi |27j. To see this, we expand Eq. (2.9) in components. One finds 
D-term mixing, V D eD^^y.\^DY, which upon electroweak symmetry breaking generates a 
Fayet-Illiopolous (FI) term for U{l)d- Such a term triggers the breaking of the dark sector 
at the GeV scale: 

<-^-'MO.)^(l^.Vf{^){^){0) . ,2.10) 
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As discussed in the previous section, the dark sector spectrum is approximately supersym- 
metric when kinetic mixing is the only form of low-energy communication between the two 
sectors. 

When produced, the GeV scale particles can decay through the kinetic mixing to SM 
particles. The dark photon, 7^, decays directly through the kinetic mixing to pairs of SM 
leptons, . If the dark Higgs, h, is too light to decay to two dark photons, it decays at one 
loop to lepton pairs. Both of these decays are prompt on galactic scales for typical values of 
the parameters: 

7rf ^ r ~ (e^aEM"^7J"^ - 10"^^ sec, 
h^tr r ~ 47r(e^a^M^/,)"^ ~ 10~Sec, (2.11) 

where for the last step we have chosen the representative values m^^,mh = 1 GeV and 
e = 5 X 10"^ 

The decay of the lightest fermion in the dark sector has important consequences for the 
astrophysical signals of our model. If the lightest fermion mixes with the dark gaugino, it 
can always decay through the kinetic mixing to the SM photon and the gravitinc[^ 7(i 7 G. 
The lifetime is found to be: 

This decay is prompt on galactic scales for low-scale SUSY breaking, and leads to a hard 
gamma ray signature. If the lightest fermion is significantly heavier than its bosonic super- 
partner, it can also decay to its superpartner and the gravitino, 7^ —>■ '-fdG, ot h hG, with 
lifetime: 



7d-7dG li67rF2i I ml 



Id ■ 

-3_/ IGeVy < 



= 3 X lO"-' sec — ^ -^-T^r 1 - — • (2.13) 
V m^, ) 1^100 TeV; V "^7./ 

^We only consider models where gravity mediation is not the dominant source of scale generation in the 
dark sector, such that ~ F/AIp < GeV. This is the case for the general framework of low-scale gauge 
mediation. 
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Due to the phase space suppression above, the decay into the dark photon is subdominant 
when the dark sector is approximately supersymmetric as in our case. Consequently DM 
decays into 7^ and h produce hard leptons, while decays into 7^ produce hard gamma rays. 
Since the DM decays into both bosonic and fermionic states in the dark sector, we are led 
to the generic conclusion that the hard lepton signals may be correlated with hard gamma 
ray signals. These signatures are studied in detail in [5U] . 

The dark spectrum and lifetimes are constrained by the requirement that the GeV scale 
cosmology is safe. We discuss the dark sector cosmology and the resulting constraints in 
section HI 



2.6. MeV Scale: Dark Matter Splitting 

It is important for DM to evade the strong constraints on direct detection mentioned in 



section 2.2 There are three possible solutions: 

1. Very small kinetic mixing, e, between the dark sector and the SM. 

2. The DM does not directly couple to '-fd ot Z. 

3. The DM multiplets are split. 

The first solution applies when DM is charged under the light dark sector. As we can see 



from equation (2.1), the DM evades direct detection if the kinetic mixing is small enough, 
e < 10^^. Interestingly, as discussed in section |4| mixing of this size may be insufficient 
to keep the dark sector in thermal equilibrium thereby interfering with the usual WIMP 
cosmology. 

The second solution can be realized by keeping the DM neutral under both the SM and 



the light gauge group. For example, in section 3.4, we consider a U{l)y, x U{l)d model where 



DM is charged only under U{1)^, which is broken at the weak scale, while U{l)d is broken at 
the GeV scale. Kinetic equilibrium is maintained between DM and the SM through double 
kinetic mixing, as we discuss in section |4j 

The third possibility is to introduce a DM splitting 6M^ > 100 keV . Indeed in such a 
case the DM couples inelastically together with an excited state, x', to the dark gauge boson. 
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7d, or Z, suppressing direct detection [23]. This bound is all that is necessary to evade the 
current constraints, but there are two special values for the splitting that are of experimental 
significance. If the splitting is of order 100 keV, the DAMA signal [12] can be reconciled 
with the bounds from other experiments through the inelastic DM scenario (iDM) [28] (see 
however 03]). If, on the other hand, the splitting is of size 5 My, > 1 MeV, it can account 
for the anomalous production of positrons observed by the INTEGRAL satellite close to the 
Galactic Center This is the eXciting DM (XDM) proposal [29] (see however [45j). If 
there are enough DM states, both scenarios can be realized. 

Suppose first, that DM is charged under Gd- Splittings with the right parametric size for 
iDM or XDM are generated by direct couplings between DM and the light Higgses 



W^S [vnN^ + y^xx) + I/split A^x/i • (2.14) 

As discussed above, we assume that 5* interacts with the SUSY breaking sector and gets 
a weak scale VEV. A^ is a singlet and stability of DM requires A^ to be heavier than x, 
Iz/A^l > bxl- -'■^ ^'^^^ limit, we integrate out A^ and find a DM splitting of size: 

5m, = ^i;^ = 100 keV f^)' ( -^X f^^)" • (2.15) 
^ AruN V 1 / I 1 GeV / V2.5 TeV/ ^ ^ 



If X is charged under the SM, the last term in Eq. (2.14) can be replaced with a coupling to 
the SM Higgs. In that case the splitting is expected to be larger. 

There is an important caveat to the above mechanism. If S gets a weak-scale F-term, the 
X scalars receive a weak-scale splitting and the dark gauge boson couples across the splitting. 
This SUSY-breaking splitting provides another mechanism for evading the constraint from 
direct detection, but the splitting is generically too large to account for iDM or XDM. If 
we wish to include these proposals, S must receive a weak-scale VEV but should have no 
F-term to leading ordeij^ Consequently, S cannot be the NMSSM singlet. We note that 



there are more options other than Eq. (2.14) for generating an MeV size DM splitting, and 

we will employ a slightly different mechanism in appendix [K\ 

Another possibility is that the DM is charged under a non-Abelian hidden sector. In 

this case, the splittings among the DM multiplet are generated radiatively after dark sector 

^An alternative possibility is to introduce another source of SUSY breaking that lifts both x scalars above 
the fermions so that the fermions constitute DM. 
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symmetry breaking [9]. In practice, non-Abelian dark sectors are more difficult to construct 
and often require elaborate Higgs sectors In the explicit models that we study below, we 
will instead focus on the simplest possibility of a f/ (1)^ hidden sector at low energies. This is 
for illustrative purposes only, and more complicated dark sectors remain a valid possibility. 



3. Models 

In this section we use the tools described above to construct four explicit models of hidden 
sector decaying DM. There are many possible models within this framework, and these 
should be viewed only as illustrative examples. The models are roughly ordered by increasing 
complexity. We begin with a minimal U{l)d dark sector which includes all of the main ideas. 
The second model embeds U{l)d into SU{2)d at the GUT scale. The SU{2)d breaking 
generates dimension 6 DM decay. For the third model, we consider DM charged under the 
SM, and find that there is always an associated hard neutrino signal. The reader who is 
primarily interested in the new correlated signals that we propose may want to skip directly 
to this model. All four models can produce hard gammas that are correlated with the 
astrophysical leptons, but only the third model also produces hard neutrinos. Finally, we 
consider a U{1)^ x U{X)d model, where no splitting is required to evade direct detection. In 
section |4| we discuss the constraints that cosmology places on these models. In appendix |A} 
we discuss some further technical details for each model. 

In the last two models the DM is not charged under the GeV-scale dark sector. Conse- 
quently no Sommerfeld enhancement is present at all, so the astrophysical constraints are 
automatically avoided. Such models are in sharp contrast to the annihilating models of [9]. 



3.1. U{l)d: The Mimmal Model 

We begin by considering the simplest possibility, Gd = U{l)d. This model captures the main 
ideas of our framework and serves as an example for the models that follow. We assume 



a kinetic mixing between U{X)d and hypercharge, as in equation (2.9). The field content is 
listed in table [T] and the setup is illustrated in Fig. [3j All fields are assumed to have canonical 
Kahler potential at the GUT scale. In order to stress the modularity of the model, we split 
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Figure 3: The setup of our minimal model. DM is charged under the hidden sector, U{l)d, 
and decays through dimension-6 GUT scale suppressed operators into the dark sector gauge 
multiplet. U{l)d kinetically mixes with hypercharge, and this kinetic mixing has three 
important effects: (i) D-term mixing causes U{l)d to break at the GeV scale, (ii) dark 
gauge bosons decay through the kinetic mixing to leptons while decays into antiprotons are 
kinematically forbidden, and (iii) DM stays in kinetic equilibrium with the SM through the 
kinetic mixing, allowing for the usual 'WIMP Miracle' cosmology (see section |4]). 





GUT 
X X Y Y 


TeV 

Xi Xi S Ni 


GeV 
h h n 




1-1 


1-10 


1 -1 



Table 1: The matter content of the U{l)d model, where i = 1,2. We stress the modularity 
of the model by grouping the fields according to their scales. 



up the superpotential into three pieces. 
The first term leads to DM decay: 

W^decay = (MquT + X)YY + MquT^X + XXlX2 ■ 



(3.1) 



(3.2) 



Integrating out the GUT scale fields generates the second operator of equation (2.2), at 
one loop [13]. Meanwhile, it is straightforward to see from the equations of motion that no 
dimension-5 decays are generated in the superpotential. 
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The second term of equation (3.1 ) determines the DM and dark sector spectrum: 



Wbm = S {yiXiXi + 1/2X2X2) + nhh . (3.3) 

We assume that S obtains a weak-scale VEV, possibly through interactions with the SUSY 
breaking sector. The different Yukawa couplings yi^2 generate masses for xi and X2 with 
m^2 > rriy.-^. Both Xi are stable on cosmological timescales and contribute to the rehc 
density, however, DM is mostly composed of X2, whose larger mass leads to a smaller anni- 
hilation cross-section and therefore to a larger abundance. The dimension-6 decay operator 



in Eq. (2.2), leads to three body decays of X2 into Xi and dark sector gauge bosons, 7^, 
and/or gauginos, 7^. 7^ and 'jd then decay to SM leptons and photons through the channels 
described in section [275l 

At the GeV scale, this model resembles the low-energy U{1) construction of Ref. fI7\ . 
The D-term mixing, described in section 2.5, generates an effective FI term for the U{l)d, 



which triggers one of the light Higgses to get a VEV at the GeV scale. Without loss of 
generality, we take h to be the one with a non- vanishing VEV. Expanding around (h), h 



and n obtain a GeV mass through the last term of Eq. (3.3). Consequently, all fields are 
lifted, forming a GeV scale mass gap. 



The last term of Eq. (3.1 ) corresponds to two copies of the splitting mechanism described 



by Eqs. (2. 14), (2. 15) 



i=l 

Splittings are generated for both Xi, evading the constraints from direct detection. The two 
splittings are of different sizes, and we note that both iDM and XDM can be incorporated 
in this model if (5M^^ ~ 100 keV and SM^^ ~ 1 MeV, or vica versa. It would be interesting 
to conduct a more detailed study of this multi-species DM model to see if indeed the two 
scenarios can be accommodated. 

That xi and X2 are long-lived follows from an unbroken Z2 x Z^, as Mqut -^00. Xi^Xi^ 



and Ni are charged under Zg, respectively. This symmetry is broken by Eq. (3.2), resulting 
in DM decays. In appendix [A], we verify that dimension-5 decays are forbidden by a GUT 
scale symmetry. 
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GUT 
H X ^ n 


TeV 

X X ^ 


GeV 
h n' Sn 




Adj □ Adj Adj 


□ □ Adj 1 


□ Adj 1 



Table 2: The matter content for the SU{2)d U{l)d model. 

We conclude by remarking that with this field content, the absence of Landau poles below 
the GUT scale places a bound on the dark gauge coupling at the GeV scale, aa < 1/30. 

3.2. SU{2)d U{l)d: GUT Scale Symmetry Breaking 

We now consider a UV completion of the previous model, by embedding U{l)d into SU{2)d 
which is broken at the GUT scale. In the following discussion, we focus on the two new 
features of this model: (i) heavy gauge bosons generate the dimension-6 DM decay, and (ii) 
the low-energy theory contains split SU{2)d multiplets. The field content is summarized 
in table [2j Again, we assume a canonical Kahler potential and group the terms in the 
superpotential according to their role, 

W = P^decay + WqUT + W^M • (3.5) 

The first term above, triggers the GUT-scale breaking SU{2)d U{l)d, 

VTdecay = f{H) + H X^ . (3.6) 

Here H = H'^T"^ is a triplet and T" = a°'/2 are the generators of SU{2)d. In most of what 
follows we suppress color indices. We take f{H) to be a potential for H with a minimum 
at {H) = Mqut:T^. Consequently, X, which is introduced to cancel SU{2) anomalies [i6] . 
obtains a GUT-scale mass and is integrated out. 

To see the effect of the breaking, we integrate out the broken SU{2)d/U{l)d generators. 
Going to the Unitary gauge and solving for the massive vector superfields, V± = Vi =F ^^2, 
one finds an additional contribution to the Kahler potential |38] . 

5K,^ = -{v\T+^;)\i\^]T-^^) , (3.7) 
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where ipi collectively denote all fields (subject to the Unitary gauge constraint), = T^±iT'^ 
are the broken generators in the corresponding representation and 

A± = V{T+,T-}i7 = M2uT- (3.8) 

Substituting the DM states, Xa = {XijX2), X° = (Xi)X2) and light Higgs, ha = (/ii,/i2) into 
Eq. (3.7) one finds the contributions, 

- [ d^O {x\x2h\h, + x\x2h\h2) +(1^2). (3.9) 



These operators are precisely of the form of the first operator in Eq. (2.2). As in the f/(l 



model, we will require m^^ > m^^. In this case, DM is mostly composed of X2 and Eq. ( |3.9 ) 



generates 3-body decay of X2 into Xi and the lights Higgses hi and /i2. We again assume 
that there is kinetic mixing between the low-energy U{l)d and hypercharge, generated by 
integrating out fields charged under both the dark sector and SM, so that D-term mixing 
generates a GeV-scale VEV for h2, which is eaten by the gauge multiplet and decays to SM 



leptons and photons as we describe in section 2.5 



In order to minimize the low-energy field content so that it matches the U{l)d model of 
the previous section, and in order to give different masses to Xi and X2i we work with split 
SU {2)fi multiplets. We can split the triplets $ and n by coupling them to H and GUT scale 
singlets, denoted by 5"$ and s„: 

IVgut = Tr [g{H) ($<! + S^^ + n'n + s„n)] . (3.10) 

For generic 5f(if)|^ the VEV of H generates GUT scale masses for all component fields except 
for one linear combination of $3 and 5*$, which we denote by 5, and one linear combination 
of rig and s„, which we denote n. The specific linear combinations that remain light depend 
on g{H). In appendix |A| we use a discrete symmetry to prove that 5* and n remain light 
and to show that dimension-5 DM decays can be forbidden for generic superpotentials. 
The low-energy theory is dictated by the superpotential terms: 

Wbm = + S^)XX + in' + Sn)h^ . (3.11) 



■^There is in general a different polynomial of H in front of each term of Eq. (3.10), which we have 
suppressed to keep our notation compact. 
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Figure 4: A model with DM charged under the SM. DM is the neutral component of a 5 + 5 
representation of SU{5)sm ^ SU{3)c x SU{2)w x DM decays through a dimension-6 

operator into the gauge multiplet of U{l)d, which kinetically mixes with hypercharge. The 
conservation of hypercharge (at the GUT scale) implies that this decay must be accompanied 
by the associated production of a neutrino. This results in a primary neutrino spectrum that 
is correlated with the leptonic cosmic rays, and will be tested by upcoming experiments such 
as IceCube/DeepCore 



After the SU (2)^ breaking splits the multiplets, the low-energy effective superpotential is of 
the form: 

W^s = S{yi XiXi + 2/2 X2X2) + nhih2 , (3.12) 

with yi couplings of order one that depend on g{H). The projection onto the light state 
S results in couplings that are not SU{2) invariant, and the TeV scale VEV of S therefore 
generates different masses for xi and X2- As before, the DM is long-lived because as Mqut 
00, there is an unbroken x Z2 symmetry, under which Xi and Xi are separately charged for 



i = 1, 2. The third term is the same as the last term of Eq. (3.3), and the low-energy dark 
sector is thus the same as the U{l)d model. It is also straightforward to induce small DM 
splittings, in order to evade the constraints from direct detection and possibly incorporate 
iDM and XDM. This is shown in appendix |Aj 

3.3. SU{5)sM X U{l)d: SM Charged-DM and Correlated Neutrinos 

We now consider a model where DM is charged under the SSM and decays through a 
dimension-6 operator into the dark sector. In this model, DM itself is not charged un- 
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GUT 
X X Y Y 


TeV 

X X S si N 


GeV 
h h n' 


SU{5)sM 


1111 

1-1 


□ □111 




111 

1 -1 



Table 3: The matter content for our model with DM charged under SU {5)sm ^ SU{3)c x 
SU{2)w X U{1)y. 

der the GeV sector, avoiding the constraints due to Sommerfeld enhancement discussed in 
section 2.2, We take % + % to be charged under an SU{5) GUT gauge group, residing in a 
5 + 5. A schematic description of the model is shown in Fig. |4| By gauge invariance, decay 
into the dark sector must be accompanied by associated SM particle production^ If one SM 
particle is produced, it must be a neutrino or Higgs. The latter produces antiprotons, which 
are constrained by PAMELA, and thus we focus on the possibility that DM decays produce 
hard neutrinos that accompany dark sector production. The discovery of such neutrinos is 
discussed in [30] . 

An important requirement for this model is that the colored partner of DM decays faster 
than the current age of the universe. This is because there are strong constraints on stable 



colored particles, as discussed in section |2.2[ These constraints are evaded if the triplet DM 
decays through a dimension-5 operator. For this model, we assume that the canonical Kahler 
potential is supplemented by one irrelevant operator, generated at the GUT scale, 

1 



DM 



(3.13) 



MquT 

where Si is a singlet with mass: m^^ < m^^ < m^^. This mechanism can be easily arranged 
since the triplet partner is expected to be heavier than the DM, due to the RG evolution of 
their masses below the GUT scale. 

We list the field content in table [3} and again we group the superpotential terms according 
to their roles: 

W = PFdecay + W^M + W^split • (3. 14) 

*Wc thank N. Arkani-Hamcd for drawing our attention to this point. 
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The first term generates dimension-6 DM decay using the same mechanism as our U{l)d 
model of section 13.11 

W^dccay = (Mgut + X)YY + MgvtXX + Xx5f. (3.15) 



Integrating out X and Y generates the third dimension-6 decay operator of Eq. (2.2) at 
one-loop: 

^/a^xS.WJ. (3.16) 

This operator results in three-body decay, with DM decaying into one neutrino or sneutrino 
and two dark gauge bosons or gauginos, which subsequently decay to SM leptons and photons 



through the operators discussed in section 2.5 



At low energies, this model resembles the constructions above: 

Wr,M = S{xx + sl)+nhh (3.17) 
We assume that S, which may be the NMSSM singlet, gets a weak-scale VEV. This generates 



a mass for the DM and the singlet Si, which plays a role in the triplet decay of Eq. (3.13). 
As in the models above, we take U{l)d to kinetically mix with hypercharge, and the D-term 
mixing generates a GeV-scale VEV for h. With no DM splitting, this model would be ruled 
out because the DM couples strongly to the SM Z boson. This constraint is evaded by 
coupling the DM to the Higgs, which generates a small splitting: 

W,pi^ = SN^ + xHdN. (3.18) 

Here is a singlet that must be heavier than x, to ensure its stability. The resulting splitting 
is too large to account for iDM or XDM. In fact, iDM is already ruled out for this model 



by the constraints from inelastic capture in the sun, as discussed in section 2^, Finally, the 
DM is long-lived due to an unbroken Z2 at the renormalizable level, under which x, x, and 
A^ are charged. 

If the DM relic density is only determined by its SU{2)w gauge interaction, its mass is 
fixed to be: ^1.1 TeV [33]. This mass is too small to fit the FERMI excess with DM 



decays [6J. Fortunately, the second operator of the above splitting mechanism, Eq. (3.18), 
opens up a new annihilation channel into SM Higgses. This raises the DM annihilation 
cross-section, allowing for heavier masses which can fit FERMI. We discuss the DM relic 
density further in section |4j 
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Figure 5: A model with double kinetic mixing. DM, X2, is charged under U{1)^, which is 
broken by a different species, Xij the TeV scale. Decays are induced by a dimension- 
6 GUT suppressed operator into the f/(l)d gauge multiplet, which kinetically mixes with 
both hypercharge and U{1)^. This double kinetic mixing is sufficient to keep DM in kinetic 
equilibrium with the SM, preserving the usual WIMP cosmology (see section |4]). There is 
no strong constraint from direct detection because the DM does not couple directly to the 
Z or U{l)d gauge boson, and therefore no DM splitting is required. 





GUT 
X X Y Y 


TeV 

Xi Xi Si 


GeV 
h h n 


U{1), 



1-1 


1 -1 





1 -1 



Table 4: The matter content for the U{1)^ x U{l)d model, where i = 1,2. 
3.4. U{1)^ X U{l)d: No Mass SpHtUng 

We now consider a model with a U{1)^ x U{l)d hidden sector. We illustrate the basic idea 
in Fig. [5] The DM, X2, is charged under f/(l)^, which is broken at the weak scale by the 
VEV of a different species xi- It decays through a dimension-6 operator into the U{l)d 
gauge multiplet. There are two advantages to this setup. First, this model does not have a 
strong constraint from direct detection, because DM does not couple directly to the Z boson 
or 7d. Therefore, unlike the previous models, no DM splitting is required. Second, there is 



no constraint from photon or neutrino measurements, as in the model of section 3.3, since 
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DM is not charged under U{l)d- Another unique feature of this model is that 2-body decays 
dominate over 3-body decays because the decay operator contains a field, Xi, which obtains 
a weak scale VEV. 

The field content of this model is listed in table HI We assume a canonical Kahler 
potential, and we group the superpotential terms according to their role, 

= Wdecay + VrDM. (3.19) 



The first term is identical to the GUT scale interactions of the U{l)d model, Eq. (3.2) 

W^decay = (MquT + X)YY + MquT^X + XXlX2 , 



generating the second decay operator of Eq. (2.2). Integrating out bifundamentals generates 
kinetic mixing between U{l)d and U{1)^, 

-|y" d^eW-^Wd, (3.20) 

of the same size as the kinetic mixing between U {l)d and hypercharge, ~ e ~ 10^'^ — 10~^. 
A mixing of this size is small enough to keep the U {l)d mass gap at a GeV, but large enough 
to keep the f/(l)x sector in thermal equilibrium with the U{l)d sector. The latter guarantees, 
through the double kinetic mixing, that U{1)^ is in thermal equilibrium with the SM. We 
discuss the cosmology of this model in more detail in section |4} 
The terms in the superpotential relevant at low energies are: 

Wnm = S2 X2X2 + SiixiXi + ^2) + nhh , (3.21) 

where S2 receives a weak scale VEV from communicating with the SUSY breaking sector, 
giving the DM a mass. Solving for the F-term of Si, one finds VEVs for xi ^-^id Xi of 
order {S2). This breaks f/(l)x at the weak scale, and the dominant DM decay is 2-body, 



with a Xi VEV insertion resulting in the operator of Eq. (2.4). DM is long-lived because as 
^GUT — ^ 00 there is an unbroken Z2 symmetry, (X25X2) —{X2,X2)- As in the previous 
models, h gets a VEV at the GeV scale due to the D-term mixing between hypercharge and 



U{1)^. The last term of Eq. (3.21) generates a GeV scale mass gap. 
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4. Cosmology of the Dark Sector 

In this section we discuss the cosmology of the dark sector and the resulting constraints on our 
framework. We find constraints on the size of the kinetic mixing between the dark sector and 
SM, e, on the DM interactions, and on the spectrum of the GeV scale states. The cosmology 
of our model resembles the cosmology of the annihilating DM framework of Ref. [9]. For 
related discussions of the cosmology of GeV scale hidden sectors, see Refs. |2S1 El EH 113 EH] • 
Below we include several new observations and a novel emphasis on the aspects of the 
cosmology that are important for decaying DM. We begin this section by discussing the relic 
density of DM, and end with a discussion on the cosmology of light dark sector fermions, 
which can decay to observable gamma rays providing a smoking gun signature of decaying 
DM [30J. 

4.I. Thermal DM Abundance 

A model of DM must of course reproduce the observed relic density, Q^h"^ ~ 0.1. The 'WIMP 
Miracle' implies that the correct abundance is achieved if DM is in kinetic equilibrium with 
the SM when it freezes out, with a WIMP cross-section, (cr^f) ~ 3 x 10^^^ cm^ s"'^. The 
same cosmology applies for decaying DM, as mentioned in the introduction, since the decay 
rate is much longer than the age of the Universe. We now discuss how our model can satisfy 
these requirements. 

DM retains the usual thermal history by interacting with dark gauge bosons which are 
in kinetic equilibrium with the SM [291 El H^- The kinetic equilibrium is maintained by 
interacting with the SM thermal bath through the kinetic mixing, 'jd'^'su ^ 7'^sm, where 
^SM denotes any relativistic SM particle with hypercharge. This reaction remains efficient 
for temperatures in the range m^^ < T^in < (e^a|j^/7r^5fy^)Mpi, where is the number of 
relativistic degrees of freedom at temperature Tkin- For the DM to be a thermal relic with 
a WIMP cross-section, Tkin must be larger than the DM decoupling temperature, T^ec — 
m^/20. The thermal history therefore places a lower-bound on the size of the kinetic mixing: 

e > 10"^ - 10"^ (4.1) 

There is tension between this constraint, and the constraint on e from direct detection when 
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DM couples elastically to the dark photon, Eq. (2.1). One way to evade the constraint of 



Eq. (4.1) is to introduce weak scale particles charged under both the dark sector and the 
SM. Particles charged under both sectors can maintain kinetic equilibrium, but they must 
be very light, (9(100 GeV), in order to do so until T^ec- Another way to alleviate this tension 
is to introduce a DM splitting, which evades the constraint from direct detection. 

The introduction of a DM splitting can change the DM annihilation cross-section in an 
interesting way. A splitting can be generated radiatively or through Yukawa interactions, as 



we discuss in section 2^ Radiative splittings are generated after breaking non-Abelian dark 
sectors with specific matter content, however such models are significantly more complicated 
to construct [3J. A simpler alternative, when the DM is charged under U{l)d, is to couple 



it directly to the light Higgses, as in Eq. (2.14). In addition to introducing splittings, these 



interactions provide the DM a direct annihilation channel into light Higgses. This Yukawa 
annihilation rate, ay, can be parametrically related to the annihilation rate into dark gauge 



bosons, (T g , as 



ay /m^\^/(5m^\^ /0.5GeV\^/ \2 / 5m^ 



ag xnT'-TdJ V "^x ^ ^ ^ v2.5 TeV/ V 100 keV 



(4.2) 



Here Sm^ is the size of the DM splitting, Eq. (2.15). We see that the Yukawa annihilation 
channel parametrically dominates the DM relic density when m^^ < 500 MeV or when the 
splitting is sufficiently large. In this regime, the DM gauge coupling must be small in order 
for the DM to have the correct relic density. This implies that non-perturbative Sommerfeld 
enhancements to the annihilation cross-section are < (9(100). Decaying DM models in this 



regime evade the constraints from the Sommerfeld enhancements discussed in section 2^ 
and annihilating models of this type cannot achieve a large enough Sommerfeld enhancement 
to fit FERMI 0. 

A similar analysis applies when DM is charged under the SM and couples to the Z, as 



in the model of section |3.3[ A splitting is required to evade the constraints from direct 
detection, which can be introduced by coupling the DM to the SM Higgs. This opens up 
a new annihilation channel of DM into SM Higgses, which allows for a larger annihilation 



cross-section and heavier DM masses, as discussed in section 3.3 



An interesting example that has no tension between the thermal history and direct de- 



tection, and does not require a DM splitting, is our U{1)^ x U{l)d model of section 3.4 
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Here, the DM is charged under f/(l)^, which is broken at the weak scale and kinetically 
mixes with the GeV-scale dark sector, U{l)d, with mixing of order ~ 10""'^. There is no 
strong constraint on this model from direct detection because DM does not couple directly 
to the Z or the light dark photon. The light dark sector, U{l)d, stays in kinetic equilibrium 
with the SM through kinetic mixing, as discussed above. The kinetic mixing between U{1)^ 
and U{l)d keeps U{1)^, and therefore the DM, in kinetic equilibrium with U{l)d through 
the interaction 'y-^h ^ 7d/i, with h corresponding to any of the light Higgses or Higgsinos 
charged under the U{l)d- The DM is thus kept in kinetic equilibrium with the SM through 
double kinetic mixing, yielding the correct relic abundanc^ 



4.2. The Lightest Dark Sector Fermion 

The dark sector may contain light particles that are long-lived. Such fields are constrained 
by cosmology, as we discuss now. There is typically no constraint on light scalars and gauge 
bosons since both can decay through the kinetic mixing with cosmologically fast timescales. 



as in Eq. (2.11). An exception to this are stable light scalars due to an unbroken discrete 
symmetry, which we discuss below. The lightest fermion, on the other hand, must decay to 
the gravitino, if kinematically allowed, which can lead to cosmologically long lifetimes. In 
what follows, we focus on the situation where the lightest fermion mixes with the gaugino, 
and we consider separately the cases where it is heaver than, approximately degenerate 
with, or lighter than the dark gauge boson. We show that the last two scenarios require the 
lightest fermion to decay to a photon and a gravitino on sub-galactic length scales, leading 
to observable gamma ray signatures [30] . 

This regime applies when there is sizeable SUSY breaking in the dark sector > GeV. 

The fermions can annihilate into dark gauge boson pairs, with cross-section cr ~ 

5f^/(87rGeV^) which leads to an abundance ^^^h"^ — 10~^. After freezeout, the fermion 

can decay to the dark gauge boson and a gravitino, which is kinematically allowed for 

^For this model, the DM does kinetically decouples from the SM during freezeout at Ty\n = ?7t-7^ ^ m^. 
After decoupling, the DM temperature scales as T = T^/T\^im but this only modifies the relic density by an 
0(1) amount. 
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low scale SUSY breaking, v -F < 10 GeV. The dark gauge boson then decays through 



the kinetic mixing to leptons. The corresponding fermion lifetime is given by Eq. (2.13 ). 
For an abundance this small, there is no constraint from BBN for electromagnetic de- 
cays [35]. There are, on the other hand, strong constraints on electromagnetic decays 
after recombination [36], however the decay discussed above always proceeds before 
recombination and hence evades the bound. An analogous discussion applies if the 
lightest dark fermion is a Higgsino that is heavier than its scalar superpartner. We 
conclude that the dark sector is not constrained by the lightest fermion when it is 
heavier than its superpartner. 

Let us now consider the regime where the dark gaugino is approximately degenerate 
with the dark gauge boson. This is the case when the dark sector spectrum is approx- 
imately supersymmetric, for instance when D-term mixing dominates as discussed in 
section 2.5 When the temperature is above m^^, the dark sector is in kinetic equilib- 
rium with the SM and the number density of the dark bosons and fermions are of the 
same order of magnitude. When the temperature drops below m^^, the dark gauge 
bosons cannot be created from the thermal bath, and they decay instantly to SM lep- 
tons through the kinetic mixing, on a timescale much faster than the Hubble rate, as 



in Eq. (2.11). The dark gauginos, on the other hand, are long lived with abundance 
controlled by their available annihilation channels. As long as m^^^ > (9(0.85) m^^, the 
finite temperature allows dark gauginos to annihilate into a dark gauge boson pairs [22j , 



with cross-section 



- O(O.l) X ^ - 10' (a,.) j ^— j . (4.3) 

where the 0(0.1) suppression results from thermal averaging, and a-^ is the DM anni- 
hilation cross-section. 

For the above parameters, the resulting relic density is ^^^h"^ — 10~^. The gauginos 



will decay to photons and gravitinos with lifetime given by Eq. (2.12). For an abun- 
dance of this size, there is no constraint from BBN on the resulting electromagnetic 
decays (see Fig. 9 from the first reference of [35]), but the gauginos must decay before 
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recombination, r < 10^^ sec, to avoid constraints from diffuse gammas [36]. Amusingly, 
tliere is a coincidence in which the time of recombination roughly equals the amount of 
time it takes light to cross our Galaxy. As a consequence, the constraint from recom- 
bination guarantees that dark gauginos produced in our galaxy decay to observable 
gamma rays. The resulting constraint on the size of the kinetic mixing is e > 10~^ and 
fixing e ~ 5 X 10~^, the constraint on the SUSY breaking scale is < 2 x 10'' GeV. A 
possible caveat in the above argument, is that by the time the fermions decouple, the 
gauge boson are already kinetically decoupled from the thermal bath. This may alter 
the final abundance by some (order one) amount. A better understanding requires 
solving the exact Boltzmann equations, which is beyond the scope of this paper. 

m^, < 0(0.85) m^. 

Lastly, we consider the regime where the dark gaugino is significantly lighter than the 
dark gauge boson, which as in the first case requires GeV-scale SUSY breaking in the 
dark sector. If m;^^ > 0.5 m^^, a gaugino pair can annihilate into one dark gauge boson, 
and an e^e~ pair, through kinetic mixing. The resulting cross-section is suppressed by 



The abundance is ^!y^h? ~ 10^, and the BBN constraint now requires t^^ < 10^ sec. 
This constraint is rather strong and can be marginally satisfied for the parameters of 



Eq. (2.12). We see that rather large kinetic mixing and a low SUSY breaking scale are 
both necessary. Again, the gaugino decays to an observable gamma ray. Finally, we 
note that when ni;^^ < m^^/2, the gauginos must annihilate into 2e+2e^, with cross- 



section suppressed by an additional e^a-Eu relative to Eq. (4.4), ruling out models 
where the lightest fermion is lighter than m^^/2. 

To summarize our findings, we see that our model is unconstrained by the lightest fermion 
when m^^ > m^^, and that otherwise cosmological constraints imply that the dark gaugino 
must decay to gamma rays with short lifetimes compared to galactic length scales, leading 
to observable gamma ray signatures. These constraints are manifest as limits on the size of 
kinetic mixing, e, and the SUSY breaking breaking scale \/F , as we discuss above. 
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We conclude this section by noting that there may be hght particles that are completely 



stable. For example, in the model of section 3.1 , h and n are stable, which follows from their 



charges under an unbroken Z2. If h and n are heavier than the dark gauge boson, they have 
a large annihilation cross-section which is parametrically similar to the heavy gaugino case 
discussed above, thus resulting in a small relic density, Qh — 10~^. On the other hand, if 
h and n are lighter than the dark gauge boson, they will have a large abundance and the 
model is excluded. In general, light fields that are stable due to discrete symmetries must 
be heavier than, and annihilate into, the unstable and lighter dark sector fields. 

5. Discussion 

The decaying DM models proposed in this paper predict a number of signals at upcoming 
experiments. The light dark sector particles can be produced in colliders, resulting in lepton 
jets, as in the annihilating models of P, 121 E] • The dark sector can also be probed at low 
energy e'^e~ colliders and fixed target experiments [4j. These direct production experiments 
have the potential to discover the dark sector, but probably cannot tell apart decaying and 
annihilating models. On the other hand, astrophysical signals can differentiate between the 
two scenarios and provide a complementary means to probe the dark spectrum [20]. As 
we discuss above, primary photons are produced when the dark gaugino is degenerate with 
or lighter than the dark photon. This results in a hard gamma ray spectrum that can 
be discovered by HESS, AGIS, and CTA and possibly FERMI |30]. Moreover, if DM is 



charged under the SM, as in the model of section 3.3, decays produce primary neutrinos, 
resulting in a hard neutrino spectrum that can be measured at upcoming experiments such 
as IceCube/DeepCore. The situation is distinct from the annihilating models. For those, 
measurements from the GC exclude the production of primary photons and neutrinos with 
sizeable branching fractions [SI El [ISl EZl • 

We conclude with two further directions that can be explored in these models. 

• It would be interesting to construct a model that is more directly related to the SUSY 
breaking sector. We have taken DM to receive a weak scale mass by coupling it to a 
singlet. Since the DM is not required to be charged under the SM or under the dark 
sector, another interesting possibility is for the DM to reside in the SUSY breaking 
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sector, for example pseudomodulus |18] . 



The U{l)d and SU{2)d models in sections 3.1 and 3.2 respectively, include two species 



of DM xi and X2- The existence of several species has several interesting implications. 
First, there can be 'Wimponium' bound states, X1X2 and X2XI) which are cosmo- 
logically long-lived. Second, it may be possible to include both the iDM and XDM 
proposals since we have shown that both species can have MeV-sized DM splittings. 
The viability of these ideas requires further study. 
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A. The Models: Superpotentials and Charges 

In each of the models of section [3} the DM can decay through dimension-6 GUT suppressed 
operators into GeV states in the hidden sector. For these models to work, it is necessary that 
renormalizable or dimension-5 operators that allow DM decays are absent. Such dangerous 
decays can be forbidden by global discrete symmetries at the GUT scale. Such symmetries 
also forbid a GUT scale mass for the DM. In this appendix we verify that the above models 
are generic and safe, by presenting such global symmetries that forbid both dangerous DM 
decays and GUT scale masses for light fields. We also collect the full superpotentials of each 
model, for easy reference. 



32 



The superpotential of our minimal U{l)d model is given by: 

W = VTdecay + VTdM + VTsplit , 

W^decay = {Mgvt + X) YY + MgvtX X + XxiX2 , 

Wy)m = S {xiXi + X2X2) + nhh , 
2 

W^spiit = ^ (5iV2 + iV,Xi/i) . (A-1) 

4 = 1 

There are several dangerous operators that are allowed by U{l)d gauge invariance. These 
include a GUT scale mass of the form, XiXji ^ TeV scale mass for the light Higgses, Shh, 
renormalizable DM decay, X2hn, and dimension-5 DM decay operators, X2hh?. All dangerous 
operators of these types are forbidden by the Zf^ x Z4 symmetry displayed in the upper left 
of table m 

The superpotential of our SU{2)d — > model is: 

W = Wdec^y + WgUT + WbM + Wsplit , 

W^decay = f{H)+HX\ 

Wgvt = Ti [g{H) {<!?^ + S^^ + n'n + Snu' + snN + Sf^N)] , 

Wt)m = + S^)xx + {n + Sn)h\ 

IVspiit = {^ + S^)NN + N{x^ + f) + Nh\ (A-2) 



Wspht is not discussed in section 3.2, and is necessary to generate a DM splitting that evades 



the constraints from direct detection, as described in section 2.2, We also add the final two 



terms to PFgut- Once H obtains a VEV, and receive GUT scale masses while the 
charged components remain light. At low energies, Wsput takes the form 

W:^ut = SiN_N+ + N+N_) + N4xl + xl) + N^{xl + x\) + N_hl + N^hl (A-3) 



where S is the light linear combination of $3 and 5"$, as in Eq. (3.12). Expanding around 



the true minimum with {S) ~ TeV and (/i2) ~ GeV, we see that has a tadpole term 
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which induces a VEV for of order GeV^/ TeV. Consequently, (A^-) contributes to the 
mass of Xi and X2, which sphts the Xt and Xi multiplets. These sphttings allow the model 
to evade the constraints from direct detection, and to possibly incorporate the iDM and/or 
XDM proposals. In the upper right of table [5} we display a Zig symmetry which forbids 
GUT scale masses for light fields and dangerous decays for DM. 

In order to avoid a Landau pole below the GUT scale, the field content of this model 
requires that ad < 1/100. If DM annihilates only into light gauge fields, a gauge coupling of 
this size is insufficient to produce the correct DM relic density. Fortunately, PVspiit introduces 
DM annihilations into the light Higgses, which can dominate the annihilation cross-section 
and lead to the correct relic density. 

SM Charged DM: SU{5)sm x U{l)d 

The superpotential and Kahler term of our model with DM charged under the SM are given 
by: 

W = Wdec^y + WoM + WsM, 
Wdecay = {M^VT + X)YY + MqvtX X + Xx^ f , 

W^M = S{xx + N^ + si) + xHdN + nhh , 

W^M = SH^Hd + 10 f5fHd + 10}H^ + -^, 

Mgvt 

K D T7^. (A-4) 

IVl GUT 

Here Wsm denotes the usual SU{5) GUT superpotential with Majorana neutrino masses 
and the NMSSM singlet for generating the /i term. Dangerous decay operators now include 
renormalizable Yukawa couplings between DM and the SM, such as lO/xif^. Such operators 
must be forbidden, and for this reason DM cannot be a fourth flavor. In table [5} we list the 
charges under a x x Zq symmetry that forbids all dangerous decays and GUT scale 
masses, where the Z^ extends the usual R-parity to the new fields. 
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f/(l)x X U{lh 

The superpotential of our U{1)^ x f/(l)d model is given by: 



M^decay = {M^vT + X) YY + MgvtX X + XxiX2 , 
Wbm = S2X2X2 + Si{xiXi + si) +nhh. (A-5) 

This model is particularly simple since no DM splitting is required to evade the constraints 
from direct detection. There is a symmetry, listed in the lower left side of table [sj that 
forbids both dangerous DM decays and GUT scale masses for the light fields. 
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Table 5: The gauge charges and example global charges for each model. Clockwise from the 
upper left, are the C/(l)d, SU{2)a U{l)d, SU{5)sm x U{l)d, and U{1)^ x U{l)d models. 
For each model, the charges forbid renormalizable and dimension 5 dark matter decays and 
GUT scale masses for hght fields. 
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